home *** CD-ROM | disk | FTP | other *** search
/ Magnum One / Magnum One (Mid-American Digital) (Disc Manufacturing).iso / d16 / winsrc.arc / MISCRES.C < prev    next >
Text File  |  1991-06-16  |  24KB  |  808 lines
  1. /*
  2.     Resident odds and ends that don't fit anywhere else.
  3. */
  4.  
  5. #include <string.h>
  6. #include <stdlib.h>
  7. #include <stdio.h>
  8. #include <ctype.h>
  9. #include <time.h>
  10. #include <stdarg.h>
  11. #include <math.h>
  12. #ifdef __TURBOC__
  13. #include <dir.h>
  14. #endif
  15. #include "fractint.h"
  16. #include "fractype.h"
  17. #include "helpdefs.h"
  18.  
  19. /* routines in this module    */
  20.  
  21. extern    void restore_active_ovly(void );
  22. extern    void findpath(char *filename,char *fullpathname);
  23. extern    void notdiskmsg(void );
  24. extern    int cvtcentermag(double *Xctr,double *Yctr,double *Magnification);
  25. extern    void updatesavename(char *name);
  26. extern    int check_writefile(char *name,char *ext);
  27. extern    int check_key(void );
  28. extern    int timer(int timertype,int(*subrtn)(),...);
  29. extern    void showtrig(char *buf);
  30. extern    int set_trig_array(int k,char *name);
  31. extern    void set_trig_pointers(int which);
  32. extern    int tab_display(void );
  33. extern    int endswithslash(char *fl);
  34. extern    int ifsload(void);
  35. extern    int find_file_item(char *filename,char *itemname,FILE **infile);
  36. extern    int file_gets(char *buf,int maxlen,FILE *infile);
  37. extern    void roundfloatd(double *);
  38.  
  39. static    void trigdetails(char *);
  40.  
  41. int active_ovly = -1;
  42. long timer_start,timer_interval;    /* timer(...) start & total */
  43.  
  44. extern char IFSFileName[80];
  45. extern char IFSName[40];
  46. extern float far *ifs_defn;
  47. extern int  ifs_changed;
  48. extern int  ifs_type;
  49.  
  50. extern char temp[], temp1[256];   /* temporary strings          */
  51.  
  52. extern int  active_ovly;
  53. extern int  xdots, ydots;
  54. extern int  dotmode;
  55. extern int  show_orbit;
  56. extern int  timerflag;
  57. extern int  debugflag;
  58. extern int  maxit;
  59. extern int  fractype;
  60. extern int  got_status,curpass,totpasses,currow,curcol;
  61. extern int  fileydots;
  62. extern int  xxstart,xxstop,yystart,yystop;
  63. extern int  display3d;
  64. extern char overwrite;
  65.  
  66. /* call next when returning from resident routine and unsure whether
  67.    caller is an overlay which has been displaced */
  68. void restore_active_ovly()
  69. {
  70.    switch (active_ovly) {
  71.       case OVLY_MISCOVL:  miscovl_overlay();  break;
  72.       case OVLY_CMDFILES: cmdfiles_overlay(); break;
  73.       case OVLY_HELP:      help_overlay();     break;
  74.       case OVLY_PROMPTS:  prompts_overlay();  break;
  75.       case OVLY_LOADFILE: loadfile_overlay(); break;
  76.       case OVLY_ROTATE:   rotate_overlay();   break;
  77.       case OVLY_PRINTER:  printer_overlay();  break;
  78.       case OVLY_LINE3D:   line3d_overlay();   break;
  79.       case OVLY_ENCODER:  encoder_overlay();  break;
  80.       case OVLY_CALCFRAC: calcfrac_overlay(); break;
  81.       case OVLY_INTRO:      intro_overlay();    break;
  82.       }
  83. }
  84.  
  85.  
  86. void findpath(char *filename, char *fullpathname) /* return full pathnames */
  87. {
  88.    if (filename[0] == '\\'
  89.      || (filename[0] && filename[1] == ':')) {
  90.       strcpy(fullpathname,filename);
  91.       return;
  92.       }
  93.    fullpathname[0] = 0;             /* indicate none found */
  94. #ifdef __TURBOC__                /* look for the file */
  95.    strcpy(fullpathname,searchpath(filename));
  96. #else
  97.    _searchenv(filename,"PATH",fullpathname);
  98. #endif
  99.    if (fullpathname[0] != 0)            /* found it! */
  100.       if (strncmp(&fullpathname[2],"\\\\",2) == 0) /* stupid klooge! */
  101.      strcpy(&fullpathname[3],filename);
  102. }
  103.  
  104.  
  105. void notdiskmsg()
  106. {
  107. static char far sorrymsg[]={"\
  108. I'm sorry, but because of its random-screen-access algorithms, this\n\
  109. type cannot be created using a real-disk based 'video' mode."};
  110.    stopmsg(1,sorrymsg);
  111. }
  112.  
  113.  
  114. /* convert corners to center/mag */
  115. int cvtcentermag(double *Xctr, double *Yctr, double *Magnification)
  116. {
  117.    extern double xxmax,xxmin,yymax,yymin,xx3rd,yy3rd;
  118.    double Width, Height, Radius, Ratio;
  119.    Width  = xxmax - xxmin;
  120.    Height = yymax - yymin;
  121.    Ratio = Height / Width;
  122.    if(xx3rd != xxmin || yy3rd != yymin || Width < 0
  123.      || (Width > 1e-8 && (Ratio <= 0.74 || Ratio >= 0.76))
  124.      || Ratio < 0.66 || Ratio > 0.84)
  125.       return(0);
  126.    /* calculate center and magnification */
  127.    Radius = Height / 2.0;
  128.    *Xctr = xxmin + (Width / 2.0);
  129.    *Yctr = yymin + Radius;
  130.    *Magnification = 1.0 / Radius;
  131.    return(1);
  132. }
  133.  
  134.  
  135. void updatesavename(char *filename) /* go to the next file name */
  136. {
  137.    char *save, *hold;
  138.    char name[80],suffix[80];
  139.    char *dotptr;
  140.  
  141.    strcpy(name,filename);
  142.    suffix[0] = 0;
  143.    if ((dotptr = strrchr(name,'.')) != NULL
  144.      && dotptr > strrchr(name,'\\')) {
  145.       strcpy(suffix,dotptr);
  146.       *dotptr = 0;
  147.       }
  148.  
  149.    hold = name + strlen(name) - 1; /* start at the end */
  150.    while(hold >= name && (*hold == ' ' || isdigit(*hold))) /* skip backwards */
  151.       hold--;
  152.    hold++;            /* recover first digit */
  153.    while (*hold == '0')         /* skip leading zeros */
  154.       hold++;
  155.    save = hold;
  156.    while (*save) {        /* check for all nines */
  157.       if (*save != '9')
  158.      break;
  159.       save++;
  160.       }
  161.    if (!*save)            /* if the whole thing is nines then back */
  162.       save = hold - 1;        /* up one place. Note that this will eat */
  163.                 /* your last letter if you go to far.     */
  164.    else
  165.       save = hold;
  166.    sprintf(save,"%d",atoi(hold)+1); /* increment the number */
  167.    strcpy(filename,name);
  168.    strcat(filename,suffix);
  169. }
  170.  
  171. int check_writefile(char *name,char *ext)
  172. {
  173.  /* after v16 release, change encoder.c to also use this routine */
  174.    char openfile[80];
  175.    char opentype[20];
  176.    int i;
  177. nextname:
  178.    strcpy(openfile,name);
  179.    strcpy(opentype,ext);
  180.    for (i = 0; i < strlen(openfile); i++)
  181.       if (openfile[i] == '.') {
  182.      strcpy(opentype,&openfile[i]);
  183.      openfile[i] = 0;
  184.      }
  185.    strcat(openfile,opentype);
  186.    if (access(openfile,0) != 0) /* file doesn't exist */
  187.       return 0;
  188.    /* file already exists */
  189.    if (overwrite == 0) {
  190.       updatesavename(name);
  191.       goto nextname;
  192.       }
  193.    return 1;
  194. }
  195.  
  196.  
  197. int check_key()
  198. {
  199.    int key;
  200.    if((key = keypressed()) != 0) {
  201.       if(key != 'o' && key != 'O')
  202.      return(-1);
  203.       getakey();
  204.       if (dotmode != 11)
  205.      show_orbit = 1 - show_orbit;
  206.    }
  207.    return(0);
  208. }
  209.  
  210.  
  211. /* timer function:
  212.      timer(0,(*fractal)())        fractal engine
  213.      timer(1,NULL,int width)        decoder
  214.      timer(2)                encoder
  215.   */
  216. int timer(int timertype,int(*subrtn)(),...)
  217. {
  218.    va_list arg_marker;    /* variable arg list */
  219.    char *timestring;
  220.    time_t ltime;
  221.    FILE *fp;
  222.    int out;
  223.    int i;
  224.    int do_bench;
  225.  
  226.    va_start(arg_marker,subrtn);
  227.    do_bench = timerflag; /* record time? */
  228.    if (timertype == 2)     /* encoder, record time only if debug=200 */
  229.       do_bench = (debugflag == 200);
  230.    if(do_bench)
  231.       fp=fopen("bench","a");
  232.    timer_start = clock();
  233.    switch(timertype) {
  234.       case 0:
  235.      out = (*subrtn)();
  236.      break;
  237.       case 1:
  238.      i = va_arg(arg_marker,int);
  239.      out = decoder(i);         /* not indirect, safer with overlays */
  240.      break;
  241.       case 2:
  242.      out = encoder();         /* not indirect, safer with overlays */
  243.      break;
  244.       }
  245.    /* next assumes CLK_TCK is 10^n, n>=2 */
  246.    timer_interval = (clock() - timer_start) / (CLK_TCK/100);
  247.  
  248.    if(do_bench) {
  249.       time(<ime);
  250.       timestring = ctime(<ime);
  251.       timestring[24] = 0; /*clobber newline in time string */
  252.       switch(timertype) {
  253.      case 1:
  254.         fprintf(fp,"decode ");
  255.         break;
  256.      case 2:
  257.         fprintf(fp,"encode ");
  258.         break;
  259.      }
  260.       fprintf(fp,"%s type=%s resolution = %dx%d maxiter=%d",
  261.       timestring,
  262.       curfractalspecific->name,
  263.       xdots,
  264.       ydots,
  265.       maxit);
  266.       fprintf(fp," time= %ld.%02ld secs\n",timer_interval/100,timer_interval%100);
  267.       if(fp != NULL)
  268.      fclose(fp);
  269.       }
  270.    return(out);
  271. }
  272.  
  273.  
  274. extern void lStkSin(void) ,dStkSin(void) ,mStkSin(void) ;
  275. extern void lStkCos(void) ,dStkCos(void) ,mStkCos(void) ;
  276. extern void lStkSinh(void),dStkSinh(void),mStkSinh(void);
  277. extern void lStkCosh(void),dStkCosh(void),mStkCosh(void);
  278. extern void lStkExp(void) ,dStkExp(void) ,mStkExp(void) ;
  279. extern void lStkLog(void) ,dStkLog(void) ,mStkLog(void) ;
  280. extern void lStkSqr(void) ,dStkSqr(void) ,mStkSqr(void) ;
  281. extern void lStkRecip(void) ,dStkRecip(void) ,mStkRecip(void) ;
  282. extern void StkIdent(void);
  283. extern void lStkTan(void) ,dStkTan(void) ,mStkTan(void) ;
  284. extern void lStkTanh(void),dStkTanh(void),mStkTanh(void);
  285. extern void lStkCoTan(void),dStkCoTan(void),mStkCoTan(void);
  286. extern void lStkCoTanh(void),dStkCoTanh(void),mStkCoTanh(void);
  287. extern void lStkCosXX(void) ,dStkCosXX(void) ,mStkCosXX(void) ;
  288.  
  289. unsigned char trigndx[] = {SIN,SQR,SINH,COSH};
  290. void (*ltrig0)() = lStkSin;
  291. void (*ltrig1)() = lStkSqr;
  292. void (*ltrig2)() = lStkSinh;
  293. void (*ltrig3)() = lStkCosh;
  294. void (*dtrig0)() = dStkSin;
  295. void (*dtrig1)() = dStkSqr;
  296. void (*dtrig2)() = dStkSinh;
  297. void (*dtrig3)() = dStkCosh;
  298. void (*mtrig0)() = mStkSin;
  299. void (*mtrig1)() = mStkSqr;
  300. void (*mtrig2)() = mStkSinh;
  301. void (*mtrig3)() = mStkCosh;
  302.  
  303. struct trig_funct_lst trigfn[] =
  304. /* changing the order of these alters meaning of *.fra file */
  305. /* maximum 6 characters in function names or recheck all related code */
  306. {
  307.    {"sin",   lStkSin,   dStkSin,   mStkSin   },
  308.    {"cosxx", lStkCosXX, dStkCosXX, mStkCosXX },
  309.    {"sinh",  lStkSinh,  dStkSinh,  mStkSinh  },
  310.    {"cosh",  lStkCosh,  dStkCosh,  mStkCosh  },
  311.    {"exp",   lStkExp,   dStkExp,   mStkExp   },
  312.    {"log",   lStkLog,   dStkLog,   mStkLog   },
  313.    {"sqr",   lStkSqr,   dStkSqr,   mStkSqr   },
  314.    {"recip", lStkRecip, dStkRecip, mStkRecip }, /* from recip on new in v16 */
  315.    {"ident", StkIdent,  StkIdent,  StkIdent  },
  316.    {"cos",   lStkCos,   dStkCos,   mStkCos   },
  317.    {"tan",   lStkTan,   dStkTan,   mStkTan   },
  318.    {"tanh",  lStkTanh,  dStkTanh,  mStkTanh  },
  319.    {"cotan", lStkCoTan, dStkCoTan, mStkCoTan },
  320.    {"cotanh",lStkCoTanh,dStkCoTanh,mStkCoTanh},
  321. };
  322. int numtrigfn = sizeof(trigfn)/sizeof(struct trig_funct_lst);
  323.  
  324. void showtrig(char *buf) /* return display form of active trig functions */
  325. {
  326.    char tmpbuf[30];
  327.    *buf = 0; /* null string if none */
  328.    trigdetails(tmpbuf);
  329.    if (tmpbuf[0])
  330.       sprintf(buf," function=%s",tmpbuf);
  331. }
  332.  
  333. static void trigdetails(char *buf)
  334. {
  335.    extern char maxfn;
  336.    int i, numfn;
  337.    char tmpbuf[20];
  338.    numfn = (curfractalspecific->flags >> 6) & 7;
  339.    if(curfractalspecific == &fractalspecific[FORMULA] ||
  340.       curfractalspecific == &fractalspecific[FFORMULA]    )
  341.       numfn = maxfn;
  342.    *buf = 0; /* null string if none */
  343.    if (numfn>0) {
  344.       strcpy(buf,trigfn[trigndx[0]].name);
  345.       i = 0;
  346.       while(++i < numfn) {
  347.      sprintf(tmpbuf,"/%s",trigfn[trigndx[i]].name);
  348.      strcat(buf,tmpbuf);
  349.      }
  350.       }
  351. }
  352.  
  353. /* set array of trig function indices according to "function=" command */
  354. int set_trig_array(int k, char *name)
  355. {
  356.    char trigname[10];
  357.    int i;
  358.    char *slash;
  359.    strncpy(trigname,name,6);
  360.    trigname[6] = 0; /* safety first */
  361.  
  362.    if ((slash = strchr(trigname,'/')))
  363.       *slash = 0;
  364.  
  365.    strlwr(trigname);
  366.  
  367.    for(i=0;i<numtrigfn;i++)
  368.    {
  369.       if(strcmp(trigname,trigfn[i].name)==0)
  370.       {
  371.      trigndx[k] = i;
  372.      set_trig_pointers(k);
  373.      break;
  374.       }
  375.    }
  376.    return(0);
  377. }
  378. void set_trig_pointers(int which)
  379. {
  380.   /* set trig variable functions to avoid array lookup time */
  381.    int i;
  382.    switch(which)
  383.    {
  384.    case 0:
  385.       ltrig0 = trigfn[trigndx[0]].lfunct;
  386.       dtrig0 = trigfn[trigndx[0]].dfunct;
  387.       mtrig0 = trigfn[trigndx[0]].mfunct;
  388.       break;
  389.    case 1:
  390.       ltrig1 = trigfn[trigndx[1]].lfunct;
  391.       dtrig1 = trigfn[trigndx[1]].dfunct;
  392.       mtrig1 = trigfn[trigndx[1]].mfunct;
  393.       break;
  394.    case 2:
  395.       ltrig2 = trigfn[trigndx[2]].lfunct;
  396.       dtrig2 = trigfn[trigndx[2]].dfunct;
  397.       mtrig2 = trigfn[trigndx[2]].mfunct;
  398.       break;
  399.    case 3:
  400.       ltrig3 = trigfn[trigndx[3]].lfunct;
  401.       dtrig3 = trigfn[trigndx[3]].dfunct;
  402.       mtrig3 = trigfn[trigndx[3]].mfunct;
  403.       break;
  404.    default: /* do 'em all */
  405.       for(i=0;i<4;i++)
  406.      set_trig_pointers(i);
  407.       break;
  408.    }
  409. }
  410.  
  411.  
  412. int tab_display()    /* display the status of the current image */
  413. {
  414.    extern char usr_floatflag;
  415.    extern double xxmin, xxmax, xx3rd, yymin, yymax, yy3rd;
  416.    extern double param[4];
  417.    extern double rqlim;
  418.    extern long calctime, timer_start;
  419.    extern int  calc_status;
  420.    extern char FormName[];
  421.    extern char LName[];
  422.    extern char IFSName[];
  423.    extern int  rseed;
  424.    extern int  invert;
  425.    int row, i;
  426.    double Xctr, Yctr, Magnification;
  427.    char msg[81];
  428.    char *msgptr;
  429.  
  430.    if (calc_status < 0)     /* no active fractal image */
  431.       return(0);        /* (no TAB on the credits screen) */
  432.    if (calc_status == 1)    /* next assumes CLK_TCK is 10^n, n>=2 */
  433.       calctime += (clock() - timer_start) / (CLK_TCK/100);
  434.    stackscreen();
  435.    helptitle();
  436.    setattr(1,0,C_GENERAL_MED,24*80); /* init rest to background */
  437.  
  438.    row = 2;
  439.    putstring(row,2,C_GENERAL_MED,"Fractal type:");
  440.    if (display3d > 0)
  441.       putstring(row,16,C_GENERAL_HI,"3D Transform");
  442.    else {
  443.       putstring(row,16,C_GENERAL_HI,
  444.        curfractalspecific->name[0] == '*' ?
  445.          &curfractalspecific->name[1] :
  446.          curfractalspecific->name);
  447.       i = 0;
  448.       if (fractype == FORMULA || fractype == FFORMULA)
  449.       {
  450.      putstring(row+1,16,C_GENERAL_HI,FormName);
  451.       i = strlen(FormName)+1;
  452.       }
  453.       trigdetails(msg);
  454.       putstring(row+1,16+i,C_GENERAL_HI,msg);
  455.       if (fractype == LSYSTEM)
  456.      putstring(row+1,16,C_GENERAL_HI,LName);
  457.       if (fractype == IFS || fractype == IFS3D)
  458.      putstring(row+1,16,C_GENERAL_HI,IFSName);
  459.       }
  460.    switch (calc_status) {
  461.       case 0:  msgptr = "Parms chgd since generated";
  462.            break;
  463.       case 1:  msgptr = "Still being generated";
  464.            break;
  465.       case 2:  msgptr = "Interrupted, resumable";
  466.            break;
  467.       case 3:  msgptr = "Interrupted, non-resumable";
  468.            break;
  469.       case 4:  msgptr = "Image completed";
  470.            break;
  471.       default: msgptr = "";
  472.       }
  473.    putstring(row,45,C_GENERAL_HI,msgptr);
  474.    if (helpmode == HELPCYCLING)
  475.       putstring(row+1,45,C_GENERAL_HI,"You are in color-cycling mode");
  476.    row += 2;
  477.  
  478.    i = 0;
  479.    if (usr_floatflag) {
  480.       putstring(row,45,C_GENERAL_HI,"Floating-point flag is activated");
  481.       i = 1;
  482.       }
  483.    if (calc_status == 1 || calc_status == 2)
  484.       if (curfractalspecific->flags&INFCALC) {
  485.      putstring(row,2,C_GENERAL_HI,"Note: this type runs forever.");
  486.      i = 1;
  487.      }
  488.    row += i;
  489.  
  490.    if (calc_status == 1 || calc_status == 2)
  491.       if (curfractalspecific->flags&NORESUME)
  492.      putstring(row++,2,C_GENERAL_HI,"Note: can't resume this type after interrupts other than <tab> and <F1>");
  493.    ++row;
  494.  
  495.    if (got_status >= 0 && (calc_status == 1 || calc_status == 2)) {
  496.       switch (got_status) {
  497.      case 0:
  498.         sprintf(msg,"%d Pass Mode",totpasses);
  499.         putstring(row,2,C_GENERAL_HI,msg);
  500.         break;
  501.      case 1:
  502.         putstring(row,2,C_GENERAL_HI,"Solid Guessing");
  503.         break;
  504.      case 2:
  505.         putstring(row,2,C_GENERAL_HI,"Boundary Tracing");
  506.         break;
  507.      case 3:
  508.         sprintf(msg,"Processing row %d (of %d) of input image",currow,fileydots);
  509.         putstring(row,2,C_GENERAL_HI,msg);
  510.         break;
  511.      }
  512.       ++row;
  513.       if (got_status != 3) {
  514.      sprintf(msg,"Working on block (y,x) [%d,%d]...[%d,%d], ",
  515.         yystart,xxstart,yystop,xxstop);
  516.      putstring(row,2,C_GENERAL_MED,msg);
  517.      if (got_status == 2) { /* btm */
  518.         putstring(-1,-1,C_GENERAL_MED,"at ");
  519.         sprintf(msg,"[%d,%d]",currow,curcol);
  520.         putstring(-1,-1,C_GENERAL_HI,msg);
  521.         }
  522.      else {
  523.         if (totpasses > 1) {
  524.            putstring(-1,-1,C_GENERAL_MED,"pass ");
  525.            sprintf(msg,"%d",curpass);
  526.            putstring(-1,-1,C_GENERAL_HI,msg);
  527.            putstring(-1,-1,C_GENERAL_MED," of ");
  528.            sprintf(msg,"%d",totpasses);
  529.            putstring(-1,-1,C_GENERAL_HI,msg);
  530.            putstring(-1,-1,C_GENERAL_MED,", ");
  531.            }
  532.         putstring(-1,-1,C_GENERAL_MED,"at row ");
  533.         sprintf(msg,"%d",currow);
  534.         putstring(-1,-1,C_GENERAL_HI,msg);
  535.         }
  536.      ++row;
  537.      }
  538.       }
  539.    putstring(row,2,C_GENERAL_MED,"Calculation time:");
  540.    sprintf(msg,"%3ld:%02ld:%02ld.%02ld", calctime/360000,
  541.       (calctime%360000)/6000, (calctime%6000)/100, calctime%100);
  542.    putstring(-1,-1,C_GENERAL_HI,msg);
  543.    row += 2;
  544.  
  545.    if (videoentry.xdots) {
  546.       sprintf(msg,"Video: %dx%dx%d %s %s",
  547.           videoentry.xdots, videoentry.ydots, videoentry.colors,
  548.           videoentry.name, videoentry.comment);
  549.       putstring(row,2,C_GENERAL_MED,msg);
  550.       }
  551.    ++row;
  552.  
  553.    putstring(row,2,C_GENERAL_MED,"Corners:                X                     Y");
  554.    putstring(++row,3,C_GENERAL_MED,"top-left");
  555.    sprintf(msg,"%20.16f  %20.16f",xxmin,yymax);
  556.    putstring(-1,17,C_GENERAL_HI,msg);
  557.    putstring(++row,3,C_GENERAL_MED,"bottom-right");
  558.    sprintf(msg,"%20.16f  %20.16f",xxmax,yymin);
  559.    putstring(-1,17,C_GENERAL_HI,msg);
  560.    adjust_corner(); /* make bottom left exact if very near exact */
  561.    if (cvtcentermag(&Xctr, &Yctr, &Magnification)) {
  562.       putstring(row+=2,2,C_GENERAL_MED,"Center: ");
  563.       sprintf(msg,"%20.16f %20.16f",Xctr,Yctr);
  564.       putstring(-1,-1,C_GENERAL_HI,msg);
  565.       putstring(-1,-1,C_GENERAL_MED,"  Mag: ");
  566.       if (Magnification < 1e6)
  567.      sprintf(msg,"%20.14f",Magnification);
  568.       else if (Magnification < 1e12)
  569.      sprintf(msg,"%20.8f",Magnification);
  570.       else
  571.      sprintf(msg,"%20.2f",Magnification);
  572.       putstring(-1,-1,C_GENERAL_HI,msg);
  573.       }
  574.    else if (xxmin != xx3rd || yymin != yy3rd) {
  575.       putstring(++row,3,C_GENERAL_MED,"bottom-left");
  576.       sprintf(msg,"%20.16f  %20.16f",xx3rd,yy3rd);
  577.       putstring(-1,17,C_GENERAL_HI,msg);
  578.       }
  579.    putstring(row+=2,2,C_GENERAL_MED,"Params,");
  580.    for (i = 0; i < 4; i++) {
  581.       sprintf(msg,"%3d: ",i+1);
  582.       putstring(-1,-1,C_GENERAL_MED,msg);
  583.       sprintf(msg,"%12.9f",param[i]);
  584.       putstring(-1,-1,C_GENERAL_HI,msg);
  585.       }
  586.  
  587.    putstring(row+=2,2,C_GENERAL_MED,"Iteration maximum: ");
  588.    sprintf(msg,"%d",maxit);
  589.    putstring(-1,-1,C_GENERAL_HI,msg);
  590.    putstring(-1,-1,C_GENERAL_MED,"     Effective bailout: ");
  591.    sprintf(msg,"%f",rqlim);
  592.    putstring(-1,-1,C_GENERAL_HI,msg);
  593.  
  594.    if (fractype == PLASMA) {
  595.       putstring(++row,2,C_GENERAL_MED,"Current 'rseed': ");
  596.       sprintf(msg,"%d",rseed);
  597.       putstring(-1,-1,C_GENERAL_HI,msg);
  598.       }
  599.  
  600.    if(invert) {
  601.       extern double f_radius,f_xcenter,f_ycenter;
  602.       putstring(++row,2,C_GENERAL_MED,"Inversion radius: ");
  603.       sprintf(msg,"%12.9f",f_radius);
  604.       putstring(-1,-1,C_GENERAL_HI,msg);
  605.       putstring(-1,-1,C_GENERAL_MED,"  xcenter: ");
  606.       sprintf(msg,"%12.9f",f_xcenter);
  607.       putstring(-1,-1,C_GENERAL_HI,msg);
  608.       putstring(-1,-1,C_GENERAL_MED,"  ycenter: ");
  609.       sprintf(msg,"%12.9f",f_ycenter);
  610.       putstring(-1,-1,C_GENERAL_HI,msg);
  611.       }
  612.  
  613.    if ((row += 2) < 23) ++row;
  614.    putstringcenter(row,0,80,C_GENERAL_LO,"...Press any key to continue...");
  615.    movecursor(25,80);
  616.    getakeynohelp();
  617.    unstackscreen();
  618.    timer_start = clock(); /* tab display was "time out" */
  619.    return(0);
  620. }
  621.  
  622.  
  623. int endswithslash(char *fl)
  624. {
  625.    int len;
  626.    len = strlen(fl);
  627.    if(len)
  628.       if(fl[--len]=='\\')
  629.      return(1);
  630.    return(0);
  631. }
  632.  
  633. /* --------------------------------------------------------------------- */
  634.  
  635. int ifsload()            /* read in IFS parameters */
  636. {
  637.    int i,j,k;
  638.    FILE *ifsfile;
  639.    char buf[201];
  640.    char *bufptr;
  641.    extern float dstack[];    /* shared temp */
  642.    int ret,rowsize;
  643.  
  644.    if (ifs_defn) { /* release prior parms */
  645.       farmemfree((char far *)ifs_defn);
  646.       ifs_defn = NULL;
  647.       }
  648.  
  649.    ifs_changed = ifs_type = 0;
  650.    rowsize = IFSPARM;
  651.    if (find_file_item(IFSFileName,IFSName,&ifsfile) < 0)
  652.       return(-1);
  653.  
  654.    fgets(buf,200,ifsfile);
  655.    strlwr(buf);
  656.    bufptr = &buf[0];
  657.    while (*bufptr) {
  658.       if (strncmp(bufptr,"(3d)",4) == 0) {
  659.      ifs_type = 1;
  660.      rowsize = IFS3DPARM;
  661.      }
  662.       ++bufptr;
  663.       }
  664.  
  665.    for (i = 0; i < (NUMIFS+1)*IFS3DPARM; ++i)
  666.       dstack[i] = 0.0;
  667.    i = ret = 0;
  668.    while (fscanf(ifsfile," %f ",&dstack[i])) {
  669.       if (++i >= NUMIFS*rowsize) {
  670.       static char far msg[]={"IFS definition has too many lines"};
  671.      stopmsg(0,msg);
  672.      ret = -1;
  673.      break;
  674.      }
  675.       }
  676.    if ((i % rowsize) != 0 || getc(ifsfile) != '}') {
  677.       static char far msg[]={"invalid IFS definition"};
  678.       stopmsg(0,msg);
  679.       ret = -1;
  680.       }
  681.    if (i == 0 && ret == 0) {
  682.       static char far msg[]={"Empty IFS definition"};
  683.       stopmsg(0,msg);
  684.       ret = -1;
  685.       }
  686.    fclose(ifsfile);
  687.  
  688.    if (ret == 0)
  689.       if ((ifs_defn = (float far *)farmemalloc(
  690.             (long)((NUMIFS+1)*IFS3DPARM*sizeof(float)))) == NULL) {
  691.      static char far msg[]={"Insufficient memory for IFS"};
  692.      stopmsg(0,msg);
  693.      ret = -1;
  694.      }
  695.       else
  696.      for (i = 0; i < (NUMIFS+1)*IFS3DPARM; ++i)
  697.         ifs_defn[i] = dstack[i];
  698.  
  699.    return(ret);
  700. }
  701.  
  702. int find_file_item(char *filename,char *itemname,FILE **infile)
  703. {
  704.    char tmpname[41];
  705.    long notepoint;
  706.    char buf[201];
  707.    int c;
  708.    if ((*infile = fopen(filename,"rt")) == NULL) {
  709.       sprintf(buf,"Can't open %s",filename);
  710.       stopmsg(0,buf);
  711.       return(-1);
  712.       }
  713.  
  714.    while (1) {
  715.       while ((c = getc(*infile)) == ' ' || c == '\t' || c == '\n') { }
  716.       if (c == EOF) break;
  717.       if (c == ';') {
  718.      while ((c = fgetc(*infile)) != '\n' && c != EOF) { }
  719.      if (c == EOF) break;
  720.      continue;
  721.      }
  722.       notepoint = ftell(*infile) - 1;
  723.       ungetc(c,*infile);
  724.       if (fscanf(*infile," %40[^ \n\t({]",tmpname) == EOF) break;
  725.       while ((c = getc(*infile)) != EOF && c != '{' && c != '\n') { }
  726.       if (c == EOF) break;
  727.       if (c == '{') {
  728.      if (stricmp(tmpname,itemname) == 0) {
  729.         fseek(*infile,notepoint,SEEK_SET);
  730.         return(0);
  731.         }
  732.      while ((c = getc(*infile)) != '}' && c != EOF) { }
  733.      if (c == EOF) break;
  734.      }
  735.       }
  736.    fclose(*infile);
  737.    sprintf(buf,"'%s' definition not found",itemname);
  738.    stopmsg(0,buf);
  739.    return(-1);
  740. }
  741.  
  742. int file_gets(char *buf,int maxlen,FILE *infile)
  743. {
  744.    int len,c;
  745.    /* similar to 'fgets', but file may be in either text or binary mode */
  746.    /* returns -1 at eof, length of string otherwise */
  747.    if (feof(infile)) return -1;
  748.    len = 0;
  749.    while (len < maxlen) {
  750.       if ((c = getc(infile)) == EOF || c == '\x1a') {
  751.      if (len) break;
  752.      return -1;
  753.      }
  754.       if (c == '\n') break;             /* linefeed is end of line */
  755.       if (c != '\r') buf[len++] = c;    /* ignore c/r */
  756.       }
  757.    buf[len] = 0;
  758.    return len;
  759. }
  760.  
  761. int matherr( struct exception *except )
  762. {
  763.     extern int debugflag;
  764.     if(debugflag)
  765.     {
  766.        static int ct = 0;
  767.        static FILE *fp=NULL;
  768.        if(fp==NULL)
  769.       fp = fopen("matherr","w");
  770.        if(ct++ < 100)
  771.        {
  772.       fprintf(fp,"err:  %d\nname: %s\narg:  %le\n",
  773.           except->type, except->name, except->arg1);
  774.       fflush(fp);
  775.        }
  776.     }
  777.     if( except->type == DOMAIN )
  778.     {
  779.     char buf[40];
  780.     sprintf(buf,"%le",except->arg1);
  781.     /* This test may be unnecessary - from my experiments if the
  782.        argument is too large or small the error is TLOSS not DOMAIN */
  783.     if(strstr(buf,"IN")||strstr(buf,"NAN"))  /* trashed arg? */
  784.                /* "IND" with MSC, "INF" with BC++ */
  785.     {
  786.        if( strcmp( except->name, "sin" ) == 0 )
  787.        {
  788.           except->retval = 0.0;
  789.           return(1);
  790.        }
  791.        else if( strcmp( except->name, "cos" ) == 0 )
  792.        {
  793.           except->retval = 1.0;
  794.           return(1);
  795.        }
  796.        }
  797.     }
  798.     return (0);
  799. }
  800.  
  801. void roundfloatd(double *x) /* make double converted from float look ok */
  802. {
  803.    char buf[30];
  804.    sprintf(buf,"%-10.7g",*x);
  805.    *x = atof(buf);
  806. }
  807.  
  808.